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ABSTRACT 

A systematic investigation of the relationship between different redshift estimation 
schemes for more than 91 000 quasars in the Sloan Digital Sky Survey (SDSS) Data Release 
6 (DR6) is presented. The publicly available SDSS quasar redshifts are shown to possess 
systematic biases of Az/(1 + z)^0.002 (600 kms^^) over both small ((5z~0.1) and large 
((5z~l) redshift intervals. Empirical relationships between redshifts based on i) Call H & 
K host galaxy absorption, ii) quasar [O ll] A3728, iii) [O ill] AA4960,5008 emission, and iv) 
cross-correlation (with a master quasar template) that includes, at increasing quasar redshift, 
the prominent Mgll A2799, Cm] A 1908 and CiV A 1549 emission lines, are estabhshed as 
a function of quasar redshift and luminosity. New redshifts in the resulting catalogue pos- 
sess systematic biases a factor of ~20 lower compared to the SDSS redshift values; system- 
atic effects are reduced to the level of Az/(1 + z)^10^'* (30kms^^) per unit redshift, or 
^2.5 X 10^^ per unit absolute magnitude. Redshift errors, including components due both to 
internal reproducibility and the intrinsic quasar-to-quasar variation among the population, are 
available for all quasars in the catalogue. The improved redshifts and their associated errors 
have wide applicability in areas such as quasar absorption outflows, quasar clustering, quasar- 
galaxy clustering and proximity-effect determinations. 
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1 INTRODUCTION 



The Sloan Digital Sky Survey (SDSS) l lYork et al.ll2000l) has pro- 
duced a revolution in both the volume and quality of spectro- 
scopic data available for quasa rs. The Data Release 5 (DR5) 
jAdelman-McCarthv et al.l 20071) and Legacy Data Release 7 
(DR7) (|Abazajian et al."2009') with their associated quasar cata- 
logues {Schneider et al. 2007, 2010, respectively) provide interme- 
diate resolution (_R~2000), moderate signal-to-noise ratio (SNR) 
(SNR~15 per 69kms~^pixel), spectra of unprecedented homo- 
geneity, covering essentially the entire "optical" wavelength region 
(A=3800-9180A). 

The quality of the Schneider et al. quasar catalogues is truly 
impressive, with errors in redshift identification reduced to the 0.01 
per cent level and individual redshift esti mates, resulting primar- 
ily from the SDSS spectroscopic pipeline jStoughton et alj[2002l . 
and the SDSS DR7 websit^, are accurate to of order A2/(l -|- 
2;)~0.002. The publication of even individual quasar redshifts, 
based on moderate resolution spectra, to such accuracy was a sig- 



nificant achievement prior to the mid-1990s, further highlighting 
the advance represented by the SDSS. 

Notwithstanding the quality of the SDSS quasar spectra 
and the associated redshift estimates, important scientific inves- 
tigations, incl uding the cl ustering of quasars themselves (e.g. 
Groom et al. 2Qol IShen e t al. 2007), the cross-correlation of 
quasars and other object pop ulations (e.g. IPadmanabhan et al] 
2009l), the proximitv effe ct (e.g. lBaitlik. Duncan. & OstrikeJI 19881 : 
Kirkman & Tytlej 20081), the origin and properties of associated 



absorbers (e.g.lNestor, Hamann. & Hidalgo"2008l: rWild et alj2003 : 
iTvtler et al .1200^ benefit significantly both from reduced systemat- 
ics in redshift determinations and the reliable assignment of redshift 
uncertainties for individual quasars. 

In this paper we present the determination of new redshifts 
and associated error estimates f or more th an 89 500 quasars from 
the SDSS DR6 ( Adelman-McC arthv et al.i r2008). Our redshift de- 
terminations suffer from much smaller systematic effects compared 
to the default values from the SDSS spectroscopic pipeline. Specif- 
ically, systematics are reduced by more than an order of magni- 
tude to 1.0x10^* in A2/(l+2) per unit redshift, or, equivalently, 
30 km s~^per unit redshifQ A detailed comparison of redshifts de- 
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The quasar research community has normally quantified redshift errors 
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rived from Call H&K absorption, [On] Willi, 31129 emission, 
[Olll] AA4960,5008 emission and cross-correlation with a new 
quasar template spectrum, provides greatly improved error esti- 
mates for individual quasar redshifts. The error estimates incor- 
porate both the uncertainties resulting from the properties of the 
SDSS spectra, quantified using the very large number of multiple 
spectra present in the SDSS, and the intrinsic quasar-to-quasar dis- 
persion. The resulting catalogue will allow significant advances in 
many studies that rely on the determination of systemic quasar red- 
shifts with small systematics and well-determined uncertainties. 

The paper is structured as follows. Section 2 describes the 
quasar sample, before the features of the quasar redshifts avail- 
able from the SDSS spectroscopic pipeline are illustrated in Sec- 
tion 3. Section 4 includes a description of the procedures involved 
in generating a master quasar template for the cross-correlation red- 
shift estimates. Section 5 then describes the procedures employed 
to provide the new redshift estimates for the SDSS quasars. An as- 
sessment of the consistency of the different redshift estimates is 
given in Section 6 and redshift estimates, based on different rest- 
frame wavelength regions, are placed onto the same 'systemic' ref- 
erence system. A critical assessment of the internal and external 
reliability of the new quasar redshifts is presented at this point. 
Twenty-one centimetre radio observations of the majority of the 
SDSS quasars are available from the Faint Images of the Radio Sk y 
at Twenty centimetres (FIRST. IsTcker, White, & Helfandl IT995h . 
Section 7 contains a description showing how the new redshift es- 
timation scheme allows spectral energy distribution (SED) depen- 
dent composite spectra (for FIRST-detected quasars in this case) 
to be constructed, producing significantly improved redshifts. The 
resulting redshift catalogue, including well-determined error esti- 
mates for each quasar, is described in Section 8. A short discus- 
sion, including consideration of the origin of the differences with 
published redshifts and an independent test of the new redshifts 
follows in Section 9. The paper concludes with a brief summary 
of the conclusions as Section 10. We adopt the same convention 
as employed in the SDSS and use vacuum wavelengths throughout 
the paper. Absolute magnitudes are calculated in a cosmology with 
Ho=70kms'\ Q.M=Q.'i and Q.a=Q.1. 



2 QUASAR SAMPLE 

The quasar sa mple consists o f 91 66 5 objects, including 77 392 
quasars in the ISchneider et alj l l2007h DR 5 catalogue that are re - 
tained in the later DR7 quasar catalogue o f ISchneider eTal] ( |201(]|) . 
A further 13 081 objects are quasars, present in the additional DR6 
spectroscopic plates, identified by on e of us (PCH) u s ing a sim- 
ilar prescription to that emp loved by ISchneider et alj ( |2007|) , all 
of which are present in the Ischneider et alT i 20101) "catalogue. An 
additional 1192 objects, which do not satisfy one, or both, of the 
emission line velocity w i dth or absolute magnitude criterion im- 
posed bv lSchneider et al] ( l2007h . are also included. While formally 
failing the 'quasar' definition of Schneider et al.'s DR5 and DR7 
compilations the objects are essentially all luminous active galactic 
nuclei (AGN). None of the results in the paper depend on the exact 
definition of the 'quasar' -sample used. 



The spectra were all processed th rough the sky-residual sub- 
traction scheme o f lWild & Hewet3 ( l2005i) , resulting in significantly 
improved SNR at wavelengths A>7200 A. The SNR improvement 
allows the important quasar rest-frame wavelength regions con- 
taining the Mgll A2799 and Cm] A 1908 emission lines, to con- 
tribute to the cross-correlation redshift determinations (Section |4!2l l 
to much higher redshifts than is possible using the original SDSS 
spectra. 

The SDSS DR6 contains a very large number of objects for 
which multiple spectra are available. For our quasar sample there 
are ~9000 independent pairs of spectra. The catalogue of spectrum 
pairs allows the accurate determination of redshift reproducibility 
as a function of SNR, redshift and cross-correlation amplitude and 
extensive use of the spectrum pairs is made to quantify the contri- 
bution of the SDSS spectra themselves to the quasar redshift errors. 



3 SDSS REDSHIFTS 

The SDSS spectroscopic pipeline (spectrold) incorporates a so- 
phisticated scheme^ for determining both the classification (star, 
galaxy, quasar,...) of the spectra and the redshifts of extragalactic 
obj ects. Cross-correlation redshift estimates are determined using 
the lTonrv & Davis j 19791) technique and a composite quasar tem- 
plate from jVanden Berk et al.ll200ll) . Emission lines are identified 
via a wavelet transform technique and an independent redshift esti- 
mate is derived using the observed-frame wavelength emission line 
locations and referen ce rest-frame emiss i on lin e wavelengths, the 
latter taken from the IVanden Berk et al" I j200lh composite quasar 
spectrum. The reference wavelength^ adopted from the quasar 
composite can differ from laboratory values due to the complex, 
often asymmetric , line profiles and apparent 'velocity shifts' of 
the lin e centroids ( lGaskellll98^ ; lTvtler & FanI 19921 ; [Richards et all 
I2OO2I) . 

The SDSS database and the individual FITS spectrum files 
contain extensive quantitative information on the determination and 
reliability of the different redshift estimates. However, the major- 
ity of researchers utilise the 'final'-redshift estimate z included in 
the SDSS SpecOb jAll table, the individual FITS spectrum file 
headers, or from the Schneider et al. catalogue^ 

If available, the cross-correlation redshift is adopted as the 'fi- 
nal' redshift for the spectrum. Some 88 per cent of the quasars 
possess redshifts derived from cross-correlation and more than a 
third of such objects also possess consistent emission line-based 
redshifts. A further 7 per cent of quasars, where no reliable cross- 
correlation redshift is available, possess redshifts derived from the 
emission lines. The remaining 5 per cent of quasars, including a 
large fraction of pathological objects and spectra of low SNR, have 
redshifts derived via manual inspection of the spectra. 



in terms of Az/(l+z), whereas researchers studying galaxies conventionally 
specify uncertainties in kilometres per second. The improvements possible 
in redshift determination made possible by the SDSS spectra are such that 
the 'kilometres per second' parameterisation is increasingly attractive and 
we specify the main results using both schemes. 



3.1 SDSS Princeton redshifts 

Independent spectrum classifications and redshift determinations, 
based on direct -fitting of template spectra to the data, have been 

^ http://www.sdss.org/dr7/dm/flatFiles/spSpec.html describes the 
spectrold FITS-file data model and lists the emission line wave- 
lengths. 

* In the Schneider et al. DR5 and DR7 quasar catalogues, catastrophic red- 
shift errors are virtually absent but otherwise the catalogued redshifts are 
the 'final'-redshift estimates from the pipeline reductions. 
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made at Princeton using the specBS cod^. Tiie redshift determi- 
nation, essentially via cross-correlation, differs from the implemen- 
tation employed in the s p e c t r o 1 d pipeline but the same compos- 
ite quasar template from lVanden Berk et al. was used. 



3.2 SDSS Redshift intercomparison 

Fig-Eshows a comparison of the SDSS final-redshifts and Prince- 
ton redshifts as a function of quasar redshifS The selection of the 
sub-sample of more than 70 000 spectra is conservative in that only 
spectra with high-confidence SDSS redshifts, where there is also 
no inconsistency between the cross-correlation and emission line 
redshift determinations, are used. The data in Fig. [T] should essen- 
tially represent an internal consistency check and the large differ- 
ences between redshifts, extending to ±5x10^'^, or 1500 kms^^, 
are surprising. Perhaps even more striking is the sequence of appar- 
ent discontinuities in the behaviour as a function of redshift. 

A second illustration of the extent of redshift-dependent sys- 
tematics comes from comparing the redshift derived from the lo- 
cation of the Mgll AA2796,2803 emission in each quasar spec- 
trum with the SDSS redshift. Fig. |2] presents the data for more 
than 60 000 spectra with SNR^IO Mgll emission line locations 
(from the SDSS spectroscopic pipelinqj). The rest-frame location 
of the Mg II emission line has been shown by many studies over 
the decades to be well-behaved and there is no reason to expect 
~500kms"^shifts over small redshift intervals, or, indeed, an ap- 
parent systematic 2x 10"'^ (600 kms^^) change in the location of 
the Mgll emission with increasing redshift of the quasars. The 
systematic redshift differences show similar patterns over the red- 
shift range common to both Fig.[T]and Fig.|2] Although somewhat 
more complex to interpret (Section|4!4](. the equivalent plot for the 
C III] emission. Fig. [3] also shows strong systematic effects as a 
function of redshift. The form and substantial amplitude of the sys- 
tematic and random differences in Figs. [T}[3] led to the initiation of 
the investigation presented here. 



4 MASTER QUASAR TEMPLATE CONSTRUCTION 

The generation of the high-SNR quasar template to be used to cal- 
culate cross-correlation redshifts begins with a sample of quasars 
at low redshifts that possess emission line-determined redshifts. A 
somewhat more involved procedure is then necessary to incorpo- 
rate additional quasars at higher redshifts into the master template. 
In this section the recipe for each element of the master template 
construction are outlined. 



4.1 Initial low-redshift quasar template 

The narrow forbidden emission lines of [Olll] AA4960,5008 are 
prominent in many quasar spectra with redshifts 2:<0.8 and a com- 
posite spectrum based on the combination of quasars with red- 
shifts determined via the location of [O III] emission forms the 
starting point for the construction of the master quasar template. 

^ http://spectro.princeton.edu/ 

^ All figures showing redshift differences between estimates zi and Z2 
have Az/(1 + z)={zi — Z2)/{1 + z) plotted as the y-axis. The choice 
of which estimate is used in the denominator is usually irrelevant given the 
scale of the plots. 

^ The SNR constraints applied to the use of emission lines refer to the sig- 
nificance of the emission line detection by the SDSS spectroscopic pipeline. 




SDSS Redshift 

Figure 1. Redshift differences, Az/(1 + z), between SDSS and Prince- 
ton pipeline reductions. Spectra plotted possess SDSS final-redshifts with 
high-confidence (zConf>0.9) derived via cross-correlation (zstatus=3 or 
4). Large differences between redshifts extend to amplitudes of ±5x 10~^ 
(±1500 kms^^). Particularly striking is the sequence of apparent disconti- 
nuities in the behaviour of the redshift differences as a function of redshift. 
Data for 69915 spectra ai'e included (503 lie outside the y-axis range plot- 
ted). 




SDSS Redshift 

Figure 2. Redshift differences, A2/(l + z), between SDSS redshifts and 
redshifts derived from the SDSS-determined Mg II emission line locations. 
Spectra plotted possess Mgll emission line SNR>I0. The solid red line, 
calculated using a 200 1 -point running median of the data-points, shows 
the form of the systematic trends with redshift. Systematic redshift differ- 
ences of ~500kms-ishifts over small redshift intervals are evident and, 
over a larger redshift interval, a prominent systematic trend of 2x 10~^ 
(600 km s^^) can be seen. Data for 60 190 spectra are included (2101 lie 
outside the y-axis range plotted). 



The ~19000 quasars with SDSS redshifts 2<0.85 are searched for 
the presence of [Olll] emission in a narrow wavelength interval 
(~100 A) corresponding to the predicted rest-frame location calcu- 
lated from the SDSS redshift. 

A 'continuum' is defined using a median filter of 21 pixels 
and [O III] AA4960,5008 emission then identified using a matched- 
filter detection schem e applied to a c ontinuum-subtracted version 
of each spectrum (e.g. Hewett et al. 1985). [O III] emission is often 
broad and frequently exhibits strong asymmetries ( Heckm an et al.l 
Il98lh . the small filter-scale adopted for the [O III] detection is cho- 
sen with the aim of isolating narrow, well-defined, peaks that may 
be present. The filter template consists of two Gaussian compo- 
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2 3 
SDSS Redshift 

Figure 3. Redshift differences, A^/(l + z), between SDSS redsliifts and 
redshifts derived from tlie SDSS-determined C III] emission line loca- 
tions. Spectra plotted possess Cm] emission line SNR>10. The wave- 
length adopted for the C III] emission centroid has been chosen to produce, 
on average, zero-offset for 2<2.2. The solid red line, calculated using a 
2001 -point running median of the data-points, shows the form of the sys- 
tematic trends with redshift. Large systematic changes of up to 3 X 10~^ 
(900kms~^) over small redshift intervals can be seen. The prominent 
horizontal 'feature' close to ~4xl0^'^ is an artifact resulting from the 
±1500 km s~^velocity interval (about the cross-correlation redshift) within 
which the SDSS reduction pipeline searches for emission lines. Data for 
59 260 spectra are included (2628 lie outside the y-axis range plotted). 

nents of the same width, centred at 4960.30 A and 5008.24 A, with 
flux ratio 1:3. 

Emission features can be identified reliably via detections with 
relatively low SNR, particularly given the restricted wavelength 
range searched in each spectrum. However, given the importance of 
establishing accurate redshifts as the first step in the construction of 
the composite quasar, the 8542 spectra possessing [O III] detections 
with SNR^8(T form the starting point for the template construction. 

The recipe used to combine spectra with specified redshifts 
into a composite is as follows: 

• pixels falling within 6.0 A of the strong night-sky lines at 
5578.5 A and 6301.7 A are flagged 

• pixels without valid SDSS data, determined from the SDSS 
noise array provided for each spectrum, are flagged 

• spectra are shifted to the rest-frame, with the native 
69 km s~^ 'pixels' of the original SDSS spectra retained. The sig- 
nal from each spectrum is placed onto the master rest-frame wave- 
length array using a 'nearest pixel'-scheme, thereby avoiding the 
need for any rebinning or interpolatior[f| 

• spectra are normalised using a wavelength interval common to 
all spectra 

• spectra are median-filtered with a window of 61 -pixels to de- 
fine a 'continuum'. Spectrum pixels falling more than 4.5(7 below 
the continuum, along with a grow-radius of two pixels, are flagged, 
effectively removing wavelengths affected by strong narrow ab- 
sorption 

• the median value of all non-flagged pixels at each rest-frame 



° The additional 'jitter' that the simple nearest-pixel scheme introduces is 
small, with a maximum error of 34.5 km s~^ and an increased dispersion of 
tT=20kms~^, less than a third of a pixel, in the extent of features in the 
resulting composite spectra. 



wavelength is calculated (a minimum of 100 spectra must con- 
tribute) 

At this point a very high SNR composite quasar spectrum 
is available, extending down to rest-frame wavelength A~2300 A. 
The [O III] emission moves beyond the red limit of the SDSS spec- 
tra at 2>0.8 and it is necessary to use a cross-correlation scheme, 
employing a much greater wavelength range of the quasar spec- 
trum, to allow the construction of the master template further into 
the rest-frame ultra-violet. 

4.2 Cross-correlation redshift algorithm 

The cross-correlation algorithm is based on a straightforward spa- 
tial cross-correlation between an individual quasar spectrum and 
a high-SNR template spectrum. The key elements of the cross- 
correlation calculation are i) a conservative choice of the portions 
of the quasar spectrum to employ in the calculation, avoiding strong 
emission lines close to the edges of the observed spectrum, and ii) 
application of an essentially identical 'window' to both the indi- 
vidual quasar and template spectra prior to the cross-correlation 
calculation. 

For each quasar spectrum, with its companion error-array, pix- 
els are excluded from the cross-correlation calculation according to 
a sequence of rules/tests. The first and last valid pixels, where the 
SDSS spectrum error-array is not set to 0, define the limits of the 
accessible wavelength range. In the observed-frame: 

• the first 25-pixels at each end are excluded 

• pixels within 6 A of each of the strong night-sky emission lines 
at 5578.5 A and 6301.5 A are excluded 

• narrow absorption features are identified by examining a 
continuum-subtracted spectrum. The continuum is defined using a 
61 -pixel median filter. Pixels that fall more than 4.5(7 below the 
continuum are flagged and a grow-radius of 2-pixels then applied. 
Thus, a single pixel exceeding the threshold results in the exclusion 
of 5-pixels. 

The quasar spectrum is then transformed to the rest-frame us- 
ing the specified redshift estimate, 2:i„ij3. In the rest-frame: 

• pixels with A>7000 A are excluded 

• for Zi„it>0.38, pixels with A>6400 A, i.e. the Hof region, are 
excluded 

• for Zinit<0A5, pixels with A<2900 A, i.e. the Mgll region, 
are excluded 

• for 2i„it<1.10, pixels with A<1975A, i.e. the Cm] region, 
are excluded 

• for Zinit <4.00, pixels with A< 1675 A, i.e. the C IV region, are 
excluded 

• pixels with A< 1275 A, i.e. the N V and Lyman-o lines and the 
Lyman-a forest, are always excluded. 

Following the definition of the restricted wavelength interval 
over which the quasar spectrum is retained, continua, estimated us- 
ing a large-scale, 601-pixel, median filter are subtracted from the 
quasar and the template spectra. Exactly the same wavelength in- 
terval is used to estimate the continua subtracted from the quasar 
and the template spectra. 

^ The SDSS-derived redshift is used to determine the value of Zi^u ini- 
tially but all the cross-correlation estimates are recalculated using an up- 
dated value of Zinit from the cross-correlation calculation itself. The cross- 
correlation estimates converge to the 10"^ level with just one iteration. 
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With continuum-subtracted quasar, Qcs, and template, Tc. 
spectra available, the cross-correlation, for lag "1", is performed: 



cc(0 = 



(1) 



where ai is the noise, as provided in the SDSS FITS-files, and the 
'cs' sub-scripts have been omitted for clarity. 

A quadratic fit is then made to the array of cc{l) values over 
the interval l=±npix, with npix=lOO. The fit is then refined, per- 
forming quadratic fits to narrower pixel intervals centred on the 
peak of the previous quadratic fit, with the final fit determined over 
an interval of l=±npix/5. The output consists of a redshift esti- 
mate, z fin , and a cross-correlation amplitude, cc,nax , in the range - 
l^ccmax ^ 1 , which parameterises the degree of similarity between 
the two spectra. Extensive experimentation demonstrates that the 
requirement cCmax^O-2 results in an almost error-free catalogue 
of cross-correlation redshifts. However, it must be stressed that the 
use of such a low amplitude is only possible because of the ex- 
tremely low occurrence of catastrophic redshift mis-identifications 
resulting from the SDSS spectroscopic pipeline and the subsequent 
refinements of .Schneider et al.i t2007i . i2010, ). 
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Figure 5. The observed-frame ratio of the centroid wavelengths for 
Mg II and C III] emission lines as a function of the quasar absolute mag- 
nitude. Mi . The quasar sample consists of 42 609 objects with redshifts 
1.1<2<2.2. The solid red line, calculated using a 2001-point running me- 
dian of the data-points, shows the strong systematic trend with quasar lumi- 
nosity. Data for 1000 spectra lie outside the x-y range plotted. 



4.3 Quasar template extension for redshifts 0.8<2:^1.6 

With the cross-correlation redshift determination procedure in 
place it is possible to utilise quasars with redshift 2>0.8 to ex- 
tend the master template. To ensure that the template construction 
is not adversely affected by the inclusion of spectra with poor SNR, 
or the presence of broad absorption line (BAL) troughs, the full 
quasar sample (Section^ was restricted to those objects satisfying 
the following criteria: 

• SDSS spectrum spectroscopic SNR, SN_ R-i-SN_I ^18. 

• quasar not identified as BAL-quasars bv lGibsonet"ai] ( l2009l) 
or from our own BAL catalogue (Allen et al. 2010, in preparation) 

Application of the criteria reduce the sample by approximately 
a half, to ~44 500 spectra. Cross-correlation redshifts are then cal- 
culated for 4071 spectra with 0.8<2:^1.0 according to the pre- 
scription of Section [431 All spectra with ccmax^O.2 and red- 
shifts, 0.8<z^l.O, are combined to produce a composite. Then, 
the original and new composites are combined by taking the aver- 
age, weighted by the relative number of spectra contributing at each 
wavelength. The effect is to determine redshifts for quasars using 
only the wavelength range where the initial (lower-redshift) com- 
posite is of high SNR. The procedure is then repeated for intervals 
of A2:=0.2 up to redshift z=l.6. Table [T] summarises the number 
of spectra, median absolute magnitudes and wavelength coverage 
for all of the composites used to generate the final master template 
spectrum. 

The key elements of the scheme are the use of wavelength 
regions ^1975 A for the calculation of redshifts of quasars up to 
0=1.6, thereby excluding the Cm] and CiV emission lines. The 
rest-frame ultra-violet region of interest is shown in Fig.|4] 

4.4 Luminosity-dependent emission line shifts 

Quasar luminosity-dependent systematic effects related to the rest- 
frame locations of Ca II absorption, [O II] , [O III] and IVIg II emis- 
sion are at or below the level of 30km s~^(Section |6]l. How- 
ever, the same is not true when considering the location of the 
Mg II emission line and the next prominent emission line complex 



of Cm] A 1908, Sim] A 1892 and Aim A 1857 as one moves fur- 
ther into the ultra-violet. Fig. [5] shows the ratio of the observed- 
frame centroidf"! of 'C III] ' and IVIg II for quasars in the redshift 
interval 1 . 1 <2:<2.2, where both lines are present in the SDSS spec- 
tra. The systematic trend of ~2x 10~^ in the wavelength ratio as a 
function of quasar luminosity translates directly into a systematic 
in Az/(l+z) where the Cm] emission line region contributes sig- 
nificantly to the cross-correlation signal. 

The luminosity-dependent change in the relative positions of 
the JVIg II and C III] emission lines presents a problem when con- 
sidering the generation of the master quasar template. Indeed, the 
large systematic means that care is needed when calculating red- 
shifts for a population of quasars with a range of luminosities. Con- 
sider the result of cross-correlating a quasar template with a quasar 
of high luminosity. The C III] emission in the quasar is at slightly 
smaller rest-frame wavelength than in the template and the result- 
ing redshift estimate will lie somewhere between redshift estimates 
based on the locations of the IVIg II and C III] lines alone. Perform- 
ing such cross-correlation redshift estimates for a sample of higher 
redshift quasars and then updating the template with their spectra 
will have the effect of biasing the profile/location of the IVIg II emis- 
sion line to smaller wavelengths. The effect is pernicious in that 
subsequent use of such a template to calculate redshifts for quasars 
where C III] is not even visible, will produce biased values because 
of the systematic change in the profile/location of IVIg II emission 
and other features in the template. Similar, even larger, luminosity- 
dependent systematic trends are also present in the relative loca- 
tions of the C III] emission complex and the C IV emission line. 

The existence of the systematic luminosity-dependent trends 
in the location of the C III] and C IV lines means that care must be 
taken in the definition of the master quasar template. 



The line centroids are generated as part of the SDSS spectrold 
pipeline. 
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4.5 Quasar template extension for redshifts 1.6<2:^2.0 

As cross-correlation redshifts are calculated using only rest-frame 
wavelengths > 1975 A for quasars with redshifts z<l.6, the current 
master template is free of the luminosity dependent systematics de- 
scribed above. For redshifts it is necessary to use the wave- 
length range including the C III] emission complex to increase the 
SNR of the cross-correlation signal. However, to preserve the form 
of the composite at wavelengths >2000 A, quasars with redshifts 
1.6^z<2.0 are incorporated into the template using only wave- 
lengths <2000A. The master quasar template is extended down 
to A=1275 A in two redshift slices (TableO, with the stricture that 
the quasars in the redshift interval 1.6^2:<2.0 are not allowed to 
contribute to the template at wavelengths ^2000 A. 

The result is a final master template in which the form of the 
spectrum at A ^2000 A is appropriate to a quasar of a particular ab- 
solute mag nitude (A/i~-260 Given that the SDSS quasars pos- 
sess a significant spread in absolute magnitude, systematic trends 
in the redshift determinations using the master quasar template are 
expected but the amplitude and form of the systematic trends are 
such that reliable corrections are possible (Section O and reliable 
redshifts can be derived for objects with redshifts up to 2=4.5. 



5 NEW QUASAR REDSHIFTS 

Redshift determinations, in decreasing order of accuracy and in- 
creasing quasar redshift, can be obtained using [Olll] emis- 
sion lines, [On] emission lines, cross-correlation including the 
Mgll emission line (Mgll.cc ), cross-correlation including the 
Cm] emission line complex (Clll]_cc ) and, finally, cross- 
correlation including the C IV emission line (C IV_cc ). For the 
cross-correlation results, empirical comparisons of redshifts de- 
rived using different rest-frame wavelength regions allow conver- 
sion relations to be derived as a function of quasar absolute mag- 
nitude and redshift. The goal is to build a redshift 'ladder' for 
quasars of increasing redshift that allows the redshift estimates to 

Absolute m agnitudes. Mi, are calculated using the prescription of 
ISchneider etaP i2007.) . 



be placed on the same underlying systemic reference system. The 
sub-Sections below consider in turn each step in the ladder. 

5.1 [O II] and [O ill] narrow emission line redshifts 

Redshifts for 13291 quasars with redshifts 2:<0.84 are available via 
the detection of [O III] AA4960,5008 emission with a SNR>6.0(t 
(Section|4T](. Similarly, detection of the [O II] AA3727,3729 emis- 
sion doublet at a SNR>6.0(t provides redshift determinations for 
an additional 3844 quasars, with redshifts 2:<1.31. 

In the case of [O II] detections, a single Gaussian, centred at 
3728.60 A, is used. The [O II] doublet consists of two components 
centred at 3727.09 A and 3729.88 A. The observed component ratio 
varies from quasar to quasar but is normally in the range 0.8:1- 
0.9:1, leading to the effective wavelength of 3728.60 A adopted. 

5.2 Cross-correlation redshifts including Mg ii 

Mgll.cc redshifts are available for a further 12289 quasars with 
redshifts 2:^ 1 . 10. The minimum rest-frame wavelength involved in 
the cross-correlation is 1975 A and systematic offsets relative to 
the emission line redshifts generated in the previous sub-Section 
are not predicted or detectable. 

In the interval 1 . 1 <z<2. 1 , the Mg II_cc redshifts involve rest- 
frame wavelengths below 1800 A and the signal is increasingly 
dominated by the C III] emission as redshift increases and the 
Mgll line shifts into the far red of the SDSS spectra. Addition- 
ally, there is also the luminosity-dependent variation in the location 
of the C III] emission complex to take into account. The amplitude 
of the systematic redshift bias is small, only ~2x 10"'* too large at 
the highest redshift 2=2.1. 

Taking the Mgll.cc redshifts and the corresponding 
Mgll emission line centroid determinations from the SDSS 
pipeline allows the dependence of the cross-correlation bias on 
redshift and absolute magnitude to be quantified. Treating the 
systematic as separable in redshift and luminosity, a sample of 
~42000 quasars, in the redshift interval 1.1<2<2.1, shows that 
linear corrections to the raw cross-correlation redshifts, with 
slopes of 1.61x10^* per unit redshift and 7.2xl0~^ per unit 
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Table 1. Quasar cross-correlation template definition parameters 



Redshift Range 


Median Mi 


Number 


FD Number 


Wavelength Coverage 

(A) 


Redshift Method 


Wavelength Contribution 

(A) 


0.0-0.4 


-22.50 


3958 


3958 


2732-8004 


[OIII] 


2732-8004 


0.4-0.8 


-23.86 


4584 


4584 


2136-6550 


[OlII] 


2136-6550 


0.8-1.0 


-24.89 


4071 


4071 


1908-5099 


Mgllxc(>1975A) 


1908-5099 


1.0-1.2 


-25.38 


4762 


393 


1732^590 


Mgllxc(>1975A) 


1732^590 


1.2-1.4 


-25.77 


5118 


375 


1589^176 


Mgn.cc 01975A) 


1589^176 


1.4-1.6 


-26.15 


5087 


341 


1466-3819 


Mgll.cc(>1975A) 


1466-3819 


1.6-1.8 


-26.45 


3882 


262 


1363-3534 


Mgllxc(>1675A) 


1363-2000 


1.8-2.0 


-26.71 


2797 


169 


1275-3275 


Mgll.cc(>1675A) 


1275-2000 




1 1.5 
HW Redshift 

Figure 6. Redshift differences, A2:/(l + z), between HW redshifts and 
redshifts derived from the SDSS-determined Mg II emission line locations. 
Spectra plotted possess Mgll emission line SNR>10. The solid red line, 
calculated using a 2001-point ruiming median of the data-points, demon- 
strates the complete removal of detectable systematic trends with redshift. 
Contrast the behaviour with that shown in Fig. |2] Data for 60 120 spectra 
are included (559 lie outside the y-axis range plotted). 
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Figure 7. Redshift differences, Az / (l+z), between connected Mg Iljcc red- 
shifts and corrected C III]xc redshifts derived using the C III] emission line 
complex alone, as a function of absolute magnitude. The solid red fine is 
a 2001-point running median of the data-points, illustrating the absence of 
any detectable systematic trends as a function of quasar absolute magni- 
tude. Data for 34 891 spectra are included (1 19 lie outside the y-axis range 
plotted). 



magnitude (reducing the raw redshift estimates as redshift and 
luminosity increase), bring any residual systematic trends down 
to the ~lxlO~® level. Note that the sense and amplitude of the 
difference between the raw Mg II_cc redshifts and the Mg II emis- 
sion line centroids are entirely consistent with the existence of 
the luminosity-dependent emission line shifts and the way that 
the master quasar template spectrum is constructed (Section [4.4t . 
Corrected Mg II_cc redshifts are available for 43 728 quasars in the 
interval 1.1<2;^2.1. 



5.3 Cross-correlation redshifts including C III] 

At redshifts z>2. \ the Mg II line no longer contributes to the cross- 
correlation redshifts and the full effect of the systematic variation 
in the rest-frame locations of the Mg II and C III] emission lines 
must be taken into account. Fortunately, an empirical determina- 
tion of the systematic differences between the corrected, un-biased, 
redshifts derived above and cross-correlation redshifts using only 
the rest-wavelength region 1675<A<2650 A, termed C IIl]_cc red- 
shifts, is straightforward to make. 

The differences between the corrected Mg II_cc redshifts and 
raw C IIl]_cc redshifts, derived using a maximum rest-frame wave- 
length of A=2650A, i.e. excluding the Mgll emission line, are 
available for ~35 000 quasars with 1.1<2:<2.0. The difference as a 



function of quasar absolute magnitude is systematic and well repre- 
sented by a linear trend; a linear fit has slope 1.67x 10^^, increas- 
ing the raw redshifts for increasing bright quasars. Application of 
the correction removes any detectable systematic effects as a func- 
tion of absolute magnitude (Fig. [TJ or redshift (Fig.[8j. The same 
correction is then applied to Clll]_cc redshifts for 13859 quasars 
with z>2.1. 



5.4 Cross-correlation redshifts including C iv 

The intention throughout is to avoid the use of the rest-frame wave- 
length region including the CiV emission line, which is known 
to show large asymmetric variations in shape, and hence of the 
line centroid. However, for 3274 quasars in the redshift inter- 
val 2.1<2:<4.5, the cross-correlation signal from the rest-frame 
A> 1675 A region is too low to produce a reliable C IIl]_cc redshift. 
For these 3274 quasars a cross-correlation redshift determination 
employing the rest-frame wavelength interval A> 1275 A (C IV_cc ) 
is made. 

There is a strong luminosity-dependent bias present due to the 
systematic variation in the location of the rest-frame C IV emission 
line centroid. The situation is complicated by the presence, in a 
large number of quasars, of significant absorption blueward of the 
C IV emission centroid, which biases the line centroid to the red. To 



© 0000 RAS, MNRAS 000, 000-000 



8 P. C. Hewett and V. Wild 



n I ' ' J ' — r — 1 — r ' r 




I I 1 ■ I.I . I . I 
L2 1.4 1.6 1.8 2 

HW Redshift 



Figure 8. Redshift dift'erences,A2/{l + z), between the con'ected 
MgII_cc redshifts and corrected ClII]xc redshifts derived, as a function 
of redshift. The solid red line is a 2001 -point running median of the data- 
points showing the absence of any detectable systematic trends as a function 
of quasar redshift. Data for 35 179 spectra are included (120 lie outside the 
y-axis range plotted). 



decouple the effects on the C IV emission line of quasar luminosity 
and the presence of absorption, a sub-sample of ~25 000 quasars 
with essentially undetectable absorption blueward of the C IV emis- 
sion line centroid is definecF^. 

An empirical correction for the systematic luminosity- 
dependent redshift bias (Fig |9ll can then be made in exactly 
the same way that the C IIl]_cc redshifts were referenced to the 
un-biased system. A two-part linear fit to the absolute magni- 
tude dependence, with slope 6.67 x 10~* for Mi<-27.0 and slope 
3.90x10^* for Mi^-27.0, provides an excellent fit to the sys- 
tematic trend. The raw C IV_cc redshift estimates are increased for 
more luminous quasars, reducing systematic redshift differences to 
undetectable levels. The amplitude of the systematic bias in the raw 
C IV_cc redshifts is large, a factor of four greater than the C III] de- 
pendence, and the correction results in a substantial reduction in 
the dispersion between the C IV_cc and un-biased redshift determi- 
nations. 

Objects with significant absorption blueward of the C IV emis- 
sion line, many of which are Broad Absorption Line (BAL) 
quasars, show an extended tail of redshift deviations to high val- 
ues. A second systematic correction, as a function of the absorp- 
tion equivalent width (AEW) is then made, using the differences 
between the corrected C IIl]_cc redshifts and the absolute mag- 
nitude corrected C IV_cc values for ~38 000 quasars. The actual 
correction applied is based on the empirically-determined median 
A2;/(l + z) versus AEW relation but the amplitude of the well- 
determined correction is closely reproduced by a linear fit with 
slope -2.5x10"^ over the 0-200 range of AEW used. The addi- 
tional quasar-to-quasar dispersion in redshift at large AEW is sig- 
nificant but only 574 quasars with AEW>20 in the final catalogue 
possess C IV_cc redshifts. 

Table |2] summarises the redshift and absolute magnitude de- 
pendent corrections made to the redshifts from the different esti- 
mation schemes in the ladder 



The absorption strength is parameterised using an integrated absorp- 
tion equivalent width (AEW), calculated over the velocity range -29 000 to 
Okms~^, relative to the predicted CiV A 1549 location 
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Figure 9. Redshift dift'erences,Az/(l + z), between the connected HW 
cross-correlation redshifts and uncorrected C IV_cc redshifts derived includ- 
ing the C IV emission line, as a function of absolute magnitude. The solid 
red line is a 2001 -point running median of the data-points. Data for 34 956 
spectra are included (1 19 he outside the y-axis range plotted). 

5.5 Additional redshifts 

Redshifts for an additional 124 quasars are available, although 
one of the strict criterion described above is not satisfied, e.g. 
SNR<6.0(7 for an emission line detection, or cCmax<Q-'2- These 
redshifts are included in the catalogue but are highlighted by the 
inclusion of a special status flag. 

Finally, the emission line detection and cross-correlation 
schemes fail to provide reliable redshift estimates for 1256 quasars. 
These objects consist primarily of a mix of pathological spectra, 
including extreme BAL quasars, and spectra of very low SNR. 
The objects are included in the redshift catalogue for complete- 
ness, with reds hifts and redshift errors taken from the Schneider et 
al. catalogues jSchneider et"ai]|2007l I20T0I) and the spectrold 
pipeline for DR6. Again, the source of the redshifts is indicated in 
the catalogue via a status flag. 



6 SYSTEMIC REDSHIFTS AND REDSHIFT 
UNCERTAINTIES 

Section [5] describes the scheme adopted to obtain redshifts using 
the most reliable estimation procedure for each quasar Redshifts 
are referenced to the zero-point provided by the location of the 
[O III] AA4960,5008 emission lines. The goal is to reduce system- 
atic errors in ^z/(l + z) to the level of ^1 x 10~* (30km s~^) per 
unit redshift and ^2.5x10"^ (8kms~^) per unit absolute mag- 
nitude. In this section the question of referencing the [O III] emis- 
sion line redshifts to the systemic system defined by the quasar host 
galaxies is considered. Starting with the comparison of absorption 
line and [O III] emission redshifts, the uncertainties in redshift esti- 
mates arising from both the intrinsic quasar-to-quasar variation and 
the reproducibility of the determinations for each technique in the 
ladder are quantified. 

6.1 [O III] and host galaxy systemic redshifts 

Redshift estimates based on the detection of photospheric absorp- 
tion from stars in the spatially averaged spectrum of the quasar host 
galaxy might be expected to provide a close to 'ideal' systemic 
redshift. Given the nature of the SDSS spectra, coupled with the 
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Table 2. Systematic redshift corrections 



KeasniTi ivieinoa 


Number 


Redshift Interval 


Redshilt CoiTectioii 

(\Az/{l + z)\ per unit z) 


A/j Correction 
{\Az/(l + z)\ per unit Mi) 


[OIII] 


13291 


<0.84 






[Oil] 


3844 


<1.35 






Mg Ilxc 


12289 


<1.1 






Mgllxc 


43728 


1.1-2.1 


1.61x10"* 


7.2x10^5 


ClII]xc 


13859 


1.1^.1 




1.67x10-* 


C IV.cc 


3274 


1.5-5.5 




6.67x10"* <-27.0 


C IV_cc 








3.90x10-* ^-27.0 



large luminosity of the quasars at rest-frame optical and near-ultra- 
violet wavelengths, detection of photospheric absorption is not pos- 
sible for the majority of objects. However, a direct comparison 
of redshifts derived from photospheric Call AA3934.8,3969.6 and 
[O III] emission is possible for a sample of objects with redshifts 
z<0.4. Generating a catalogue of Call absorption detections with 
SNR>6(7, matched to quasars with [O III] detections in the redshift 
interval 0.2^2^0.4, produces 825 quasars. Restricting the sample 
to objects that satisfy the luminosity-criterion for inclusion in the 
Schneider et al. SDSS quasar compilations, results in 615 quasars 
with absolute magnitudes covering the full range -24.0<A/i<- 
22.0. Composite spectra with median A/i=-22.2 and -22.6 possess 
both Ca II absorption and [O III] emission at high SNR. 

Measuring the centroid of the strong Call K-line at 3934. SA 
(Ca 11 H is blended with He absorption, producing a shift to longer 
wavelengths) and the centroid of the [Olll] AA4960.30,5008.24 
emission, measured above the 50 per cent peak-height level, shows 
that the [Olll] emission is shifted by 45±5kms-^to the blue, 
with no detectable dependence on luminosity. The offset deter- 
mined from the composite spectra is in good agreement with the re- 
sults of ( lBorosonll2005l) and the distribution of individual Ca II and 
[Olll] redshift-differences for the 615 objects contributing to the 
composites. The median redshift difference for the sample corre- 
sponds to a velocity shift of 38±9 km s^^blueward for [O III] rela- 
tive to Ca II . We therefore correct all redshifts by a 45 km s-^shift 
to the red to bring the zero-point into coincidence with the system 
defined by the Ca II K-line. 

After correcting for the contribution of the redshift repro- 
ducibility (crrepro=2.1 X 10"*, Calculated using multiple spectra 
(Section |2}) in estimating the Call redshifts, the empirically de- 
termined quasar-to-quasar scatte^3 of [O III] redshifts about the 
Call absorber redshifts is cri„irm=3.5x 10~*. 



6.2 [O 111] and [O ii] emission line redshifts 

A similar comparison can be made between the [Olll] and 
[On] emission line derived redshifts for the more than 7500 
quasars, redshifts z<0.8, with both [Olll] and [Oil] emission 
line detections. A small systematic velocity offset is present, with 
the [On] derived redshifts 24±5kms-^redward of the [Olll] de- 
rived redshifts. Thus, relative to the systemic reference defined by 
the Call K absorption, the offsets are 21±5 kms-^blueward for 



" All root-mean-square (rms), or a, values are calculated using abso- 
lute differences, |3;i-median(a;)j, with an iterative rejection scheme that 
removes values >4a, up to a maximum of three iterations. In fact, the pa- 
rameter distributions are not in general significantly non-Gaussian, with a 
maximum of ~5 per cent of values excluded from the final rms-estimates 



[On] and 45±5kms- blueward for [Olll] , in excellent agree- 
ment with previous work (Boroson 2005, e.g.). 

Comparison of 1378 (1103) spectrum pairs results in median 
errors of 5 x 10~^ and 1.4x 10~* in the reproducibility of Az/{l+z) 
for [O III] and [O II] respectively. The smaller error for [O III] re- 
sults from the typically higher SNR of the emission compared to 
the [O II] line. 

Systematic trends, derived from linear fits to the red- 
shift differences between [Olll] and [Oil] redshifts, are 
A2/(l-l-2)=2.1xl0"'^ per magnitude and A2/(l-l-2:)=5.9x 10"^ 
per unit redshift. The sense is that redshifts derived from the 
[Olll] emission lines become systematically smaller for objects 
with higher luminosities(redshifts). The luminosity-dependent sys- 
tematic is stronger than for redshift given the dynamic ranges of the 
variables present in the sample. The trend is consistent with an in- 
creasing degree of blue-asymmetry present in the [O 111] emission 
lines at increasing quasar luminosity. The effect is, however, small, 
amounting to a maximum of ±4.5x 10~^, or, ±14kms-^, about 
the mean relation for the sample. 

The empirically determined quasar-to-quasar rms-scatter of 
the individual [Oil] redshifts about the [Olll] redshifts is 
o-i„trm=l-OxlO"* (30kms-^). 



6.3 [O III] and Mg ii_cc redshifts 

[O III] emission line redshifts compared to Mg ll_cc redshifts, 
calculated excluding the [O III] emission lines from the quasar 
template, for ~ 12 500 quasars with [Olll] emission line reshifts 
shows an undetectable offset, median | A2;/(l-l-2)|<10~^. System- 
atic trends, from linear fits to the redshift differences, in sense 
[O III] - Mg II_cc redshifts, are A2:/(l-l-2:)=-l-l x 10^^ per magnitude 
and A2/(l-l-2)=-l.l X 10~* per unit redshift. Both luminosity and 
redshift parameterisations lead to systematics of at most ±3 x 10^^ 
(lOkms"^), for the dynamic ranges present in the sample, more 
than an order of magnitude below the uncertainties in the individ- 
ual [Olll] redshifts (Section [6. lb . The sense and amplitude of the 
small systematic is consistent with the results from the [O 111] to 
[On] emission line comparison (Section [6. 2t and is again almost 
certainly due to the increasing degree of blue- asymmetry present in 
the [O III] emission lines at increasing quasar luminosity. The com- 
parison shows that the Mg ll_cc redshifts are tied to the reference 
[O III] emission line redshifts to very high accuracy and that any 
systematics present are at most 10~* in Az/(l+z), or 30kms-^in 
velocity. 

The empirically determined quasar-to-quasar scatter of 
the individual Mgll.cc redshifts about the [Olll] redshifts is 
(Ti„trin=2.5x 10~* or 75kms-^, which represents an improve- 
ment of a factor of ~3 compared to careful determinations of 
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Table 3. Redshift uncertainties 



Redshift Method 


Number 


Redshift Interval 


Reproducibility 


Population rms 


Cumulative Population rms 








{Az/(l + z)) 


(Az/(l + z)) 


(A2/(1+Z)) 


CallK 




0.2-0.4 


2.1x10-4 


- 


- 


[OIII] 


13291 


<0.84 


5x10-5 


3.5x10-* 


3.5x10-* 


[Oil] 


3844 


<1.35 


1.4x10-* 


1.0x10-* 


3.65x10-* 


Mg II.cc 


12289 


<1.1 


1.4x10-* 


2.5x10-* 


4.3x10-* 


Mg II_cc 


43728 


1.1-2.1 


3.6x10-* 


2.5x10-* 


4.3x10-* 


C III]_cc 


13859 


1.1^.1 


6.5x10-* 


2.5x10-* 


4.3x10-* 


C IV_cc 


2700 


1.5-5.5 


5.1x10-* 


8.0x10-* 


9.1x10-* 


C IV.CC + AEW 


574 






1.0x10-3 


1.4x10-3 


extra_cc 


124 


0.5-^.5 


6.0x10-* 


3.0x10-* 


3.5x10-* 


SDSS 


1256 


0.3-5.5 


as per SDSS 
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Figure 10. Redshift differences, A2/{1 + z), between the corrected 
C III]xc redshifts and corrected C IV_cc redshifts derived using the 
C IV emission line complex alone, as a function of redshift. Data for 25 008 
quasars with AEW 5520 are included (60 lie outside the y-axis range plot- 
ted). The solid red line is a 2001-point running median of the data-points. 
The offset between the two redshift estimators is minimal except for the 
systematic 'dip' centred at Z2;2.8, coincident with the significant drop in 
detection efficiency for non-BAL quasars in the SDSS. 



the M g II emission line location (e.g. iNestor. Hamann. & Hidalg3 
l2008h . 



6.4 Mg ii.cc and C iii].cc redshifts 

The luminosity-dependent correction to bring the Mgll.cc and 
Clll]_cc redshifts into coincidence is highly successful, as evi- 
denced by Fig.|7]and Fig. [8] 

After allowing for the uncertainty in the determination of the 
Mg ll.cc and C lll]_cc redshifts due to the limited SNR of the SDSS 
spectra there is no evidence for an additional quasar-to-quasar red- 
shift scatter associated with the use of the C III] emission line re- 
gion alone. 



6.5 C iii]_cc and C iv_cc redshifts 

The diversity of the form of the C IV emission line has been 
known through many studies going back decades. The removal 
of the systematic quasar luminosity-dependent behaviour, amount- 
ing to ~650kms-^(over four magnitudes in quasar absolute mag- 
nitude), improves ClV_cc redshifts considerably. However, even 
for quasars with small AEW values the empirically determined 




HW Redshift 

Figure 11. Redshift differences, Az/(1 + z), between the final HW- 
redshifts and SDSS-redshifts, as a function of redshift. Note the large range 
on the y-axis. Data for 90409 quasars are included (1355 lie outside the 
y-axis range plotted). 



quasar-to-quasar scatter of the individual C IV_cc redshifts about 
the C IIl]_cc redshifts is ai,itrin=8.0x 10^*. The situation is much 
worse for quasars with significant absorption blueward of the 
Civ emission line centroid. The systematic correction applied 
reaches a full SxlQ-^ for the most affected quasars and the dis- 
persion at fixed AEW value adds an additional rms-scatter of 
a-intrin=l X \Q~'^ . The number of quasars with large AEW values 
where only C IV_cc redshifts are available is small, just 574 objects, 
but the associated redshift uncertainty is accordingly large. 

Fig. [To] shows the redshift differences between the corrected 
Clll]_cc and corrected Civ.cc redshifts for 25 008 quasars with 
only modest absorption (AEW^20) blueward of the C IV emission 
line. The small amplitude of the systematic differences illustrates 
the success of the statistical correction to the initial C IV_cc red- 
shifts. However, the systematic negative trend centred on z~2.8, 
which reaches ~-1.8x 10^* (55 km s^^), illustrates the limitations 
of the correction. The SDSS quasar selection experiences a dra- 
matic drop in efficiency over the interval z~2.7-2.9 that is much 
greater for non-BAL quasars than for BAL quasars. Thus, the over- 
all number of quasars drops significantly, while the BAL-fraction 
increases significantly. There are only 13 quasars with 2.7^z^2.9 
and with AEW^20, in the final catalogue but a systematic of at 
least -2x10-* may well also be present among the 574 quasars 
with C IV_cc redshifts and AEW>20 
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6.6 Summary 

All redshift estimates, except those taken direct from the SDSS, 
have been increased by 45 km s^^ (Section 16. It to bring the 
[O III] emission line based estimates onto the systemic system de- 
fined by the Ca II K absorption. 

Based on the large sample of repeat spectra, the internal re- 
producibility of the new cross-correlation redshifts represent an im- 
provement of more than a factor 2 over the SDSS redshift values. 
The reproducibility of the new redshifts is indistinguishable from 
that for the Princeton redshift values up to redshift z~ 1 .6. At higher 
redshifts the Princeton algorithm employs more information, via in- 
clusion of the C IV emission line region at A< 1675 A, and results in 
significantly better reproducibility for redshifts z>2.0 than in the 
scheme presented here. However, the large quasar-to-quasar vari- 
ation contributing to the redshift uncertainty (Table [3} means that 
differences in the internal reproducibility are not a critical factor 
for the redshifts in the final catalogue. 

Table [3] summarises the uncertainties for the different redshift 
estimates. Accurate redshift reproducibility estimates are available 
for all quasars, based on emission line SNR or ccmax value, and 
are incorporated in the redshift errors included in Table |4] To pro- 
vide an indication of the relative contributions of the internal and 
quasar-to-quasar uncertainties. Col. 4 of Table [3] lists the median 
internal error for each redshift estimate. Col. 5 gives the quasar-to- 
quasar error. The quasar-to-quasar errors, working from the chosen 
Ca II absorption reference, have been added in quadrature to pro- 
duce the cumulative quasar-to-quasar uncertainties in Col. 7. 



7 QUASARS WITH DETECTIONS IN FIRST 

The procedures described in Sections |5] and [6] reduce systematic 
redshift errors as a function of redshift and absolute magnitude 
by more than order of magnitude compared to the publicly avail- 
able SDSS redshifts. However, a further significant reduction in 
the remaining relatively large quasar-to-quasar uncertainties will 
require a detailed investigation of the spectral energy distribution 
(SED) dependent changes in the properties of the most promi- 
nent emission lines (which dominate the redshift determinations). 
Such an investigation is beyond the scope of this paper but it is 
relatively straightforward to consider systematic redshift differ- 
ences that correlate with the detection of SDSS quasars in the 
Faint Images of the Radio Sk y at Twenty centimetres (FIRST, 
iBecker, White. & Helfandlll995l) . 

For redshifts z<l.l, the systematic differences between 
the populations of FIRST-detected (FD) and not-FIRST-detected 
(nFD) quasars with Mg II_cc cross-correlation redshifts are at 
the A2/(l-l-2)^10~'' level. However, for redshifts 2>1.1, once 
the C III] -l-Si III] -l-Al III emission line complex contributes to 
the redshift determination, systematic differences become increas- 
ingly evident, reaching an amplitude of nearly Az/il+z)=2x 10"'^ 
(600km s"^) at redshifts z>4. 

The origin of the redshift differences is primarily a system- 
atic change in the ratio of the C III] and Si III] emission lines in 
the FD- and nFD-populations. The line ratio change results in a 
shift in the centroid of the blended line; Si III] is weaker in the FD- 
detected quasars and the blended line centroid moves redward. The 
Mg II_cc redshifts for the FD qu asars are thus too large. 

Based on the prescription o f ISchneider etaP j2007h for match- 
ing SDSS quasars to the FIRST survey there are 4326 FD-quasars 
with 2>1.1 in the DR6 quasar catalogue. The small fraction (7 
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Figure 12. Rest-frame spectra including the prominent emission line com- 
plex of C III] A1908, Si III] A1892 and Al III A1857, for the master quasai- 
template (red dash-dot line) and the FIRST-detected template (blue solid 
line). The reference wavelengths of the three emission line species are indi- 
cated by the vertical dotted lines. The significant difference in flux coinci- 
dent with the location of the Si III] A 1892 is evident. 



per cent) of FD-quasars, combined with the low amplitude of the 
systematic redshift differences, means that the inclusion of FD- 
detected quasars in construction of the master template results in 
changes to cross-correlation redshifts of A2/(l-l-2)< 10"^. How- 
ever, generation of an individual template for the FD-quasars re- 
sults in a quasar template with significant differences in the form 
of the C III] -l-Si III] -l-Al III emission line complex (Fig. [12}. 

Construction of the new FD-quasar template proceeds in an 
identical fashion to that described in Section |4] but the individual 
quasars used differ. For redshifts 2;<1.0 the same quasars used to 
generate the master quasar template are employed, whereas, at red- 
shifts z> 1.0, only FD-quasars are used, with the minimum num- 
ber of spectra required to generate a composite in a redshift slice 
reduced from 100 to 50 (Section |4ll. The number of quasars con- 
tributing at each redshift interval is given in Col. 4 of Table [T] As 
evident from Fig.[T2l the form of the C III] -l-Si III] -l-Al III emission 
line complex differs between the master- and FD-quasar- templates. 
The size of the empirical transformations necessary to bring the 
FD-quasar redshift estimates onto the reference system, in which 
the Mg II emission line centroid does not vary, with either redshift 
or absolute magnitude, are significantly reduced compared to the 
quasar population as a whole. 

For redshifts 1.1^2<2.1 a reduction of A2/(l + 
2:)=1.42x 10""* in the Mgll.cc redshifts, independent of red- 
shift and absolute magnitude, is necessary. For C IIl]_cc redshifts 
an absolute magnitude dependent correction of 4.19x 10^^, in the 
opposite sense to that for the master template, is required, i.e. for 
bright quasars the C IIl]_cc redshifts are reduced. However, given 
the ~4mag dynamic range of the quasar sample, the maximum 
correction for any quasar is Az/(1 + 2:)<10^*. For the small 
number of quasars where it is necessary to employ the C IV emis- 
sion line, the C IV_cc redshifts require a correction with slope 
3.44xl0~*, in the same sense as the larger correction necessary 
for the master template, i.e. for bright quasars the C IV_cc redshifts 
are increased (Fig.|9]l. 
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Table 4. Quasar redshift catalogue. The full table is available in the electronic version of the journal. 



Name 


RA 

J2000 (deg) 


Dec 

J2000 (deg) 


z 




FIRST 


2 Alternate 


Plate 


MJD 


Fibre 


Code 


SDSS J000006.53+003055.2 


0.02723 


0.51534 


1.823154 


0.001025 


-1 


-999.0 


685 


52203 


467 


3 


SDSS J000008.13+001634.6 


0.03390 


0.27630 


1.836332 


0.000614 


-1 


-999.0 


685 


52203 


470 


3 


SDSS J000009.26+151754.5 


0.03860 


15.29848 


1.197436 


0.000369 





1.200035 


751 


52251 


354 


2 


SDSS J000009.38+135618.4 


0.03909 


13.93845 


2.240486 


0.001468 





2.240486 


750 


52235 


82 


7 


SDSS J000009.42-102751.9 


0.03927 


-10.46443 


1.851731 


0.000966 


-1 


-999.0 


650 


52143 


199 


3 


SDSS JOOOOll.41+145545.6 


0.04755 


14.92935 


0.460127 


0.000357 





-999.0 


750 


52235 


499 


1 


SDSS JOOOOll. 96+000225.3 


0.04984 


0.04036 


0.478321 


0.000358 


-1 


-999.0 


387 


51791 


200 


1 


SDSS J000012.25-003220.5 


0.05108 


-0.53905 


1.437047 


0.000697 


-1 


-999.0 


1091 


52902 


129 


3 


SDSS J000013.14+141034.6 


0.05479 


14.17630 


0.949947 


0.000521 





-999.0 


750 


52235 


98 


3 


SDSS J000013. 80-005446.8 


0.05751 


-0.91300 


1.840606 


0.000789 


-1 


-999.0 


1091 


52902 


108 


3 



Table 5. Quasar cross-correlation templates. The full table is available in the electronic version of the journal. 



Wavelength 


Relative Flux^ 


Number 


Relative Flux'' 


Number 


(A) 


(per unit wavelength) 




(per unit wavelength) 




1275.26 


8.130 


179 


-999.0 





1275.56 


7.582 


186 


-999.0 





1275.85 


7.845 


191 


-999.0 





1276.15 


7.780 


196 


-999.0 





1276.44 


7.673 


204 


-999.0 





1276.73 


7.682 


212 


-999.0 





1277.03 


7.730 


219 


-999.0 





1277.32 


7.862 


229 


-999.0 





1277.61 


7.644 


234 


-999.0 





1277.91 


7.755 


247 


-999.0 






Data for master template 
^ Data for FD-quasar template 

8 THE REDSHIFT CATALOGUE 

Table |4] includes redshifts and error estimates for 91 665 quasars. 
Col. 1: is the SDSS coordinate object name, taken from the SDSS 
DR7 Legacy Release whenever available. Cols. 2 and 3 give the ob- 
ject J2000 right ascension and declination in decimal degrees. The 
redshift and redshift error are given in Cols. 4 and 5 respectively. 
Col. 6 provides a code specifying the FIRST-detection (FD) status 
of the quasar (-1: not detected, 0: outside FIRST footprint, 1: de- 
tected). Col. 7, alternate redshift for quasars with Col. 6 detection 
code =0 (derived using the FD-quasar template) and =1 (derived us- 
ing the master quasar template). The alternate redshift is assigned 
a value of '-999. 0' for quasars with Col. 6 detection code =-1. 
Col. 8 specifies the origin of the redshift estimate via a numerical 
code (1: [Olll] , 2: [On] , 3: MgII_cc , 4: Clll]_cc , 5: Civ.cc , 
6: extra.cc, 7: SDSS). The SDSS spectrum from which the redshift 
estimate is derived is specified via the spectroscopic plate number, 
modified Julian date of observation and fibre number in Cols. 9, 10 
and 1 1 respectively. The redshifts are given to six decimal places 
but, as evident from the size of the associated errors, the accuracy 
for individual objects is two orders of magnitude larger. The high 
level of precision is retained to avoid quantisation when comparing 
different redshift estimates specified to only four decimal places. 

The provision of alternate redshifts for FD-quasars and 
quasars whose FIRST detection status is unclear, allows the use 
of an appropriate redshift by researchers with particular definitions 
of 'radio' -quasar subsamples and/or additional radio-observations 
for quasars outside the current FIRST footprint. The primary (Col. 
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4) and alternate (Col. 7) redshifts differ only when the primary red- 
shift is derived from cross-correlation (with one of the two quasar 
templates) and has a value z>\.\. 

Two quasars, SDSS J134415.75-^331719.1 and 
SDSSJ142507.32-I-323137.4, exhibit distinctive double-peaked 
narrow emission. In both cases, the redshift corresponding to the 
higher velocity system is included in the table. 

The majority of researchers will be interested in the combined 
redshift error (Col. 5) arising from the limited SNR of the SDSS 
spectra and intrinsic variation from quasar-to-quasar. However, the 
internal contribution can be recovered straightforwardly via use of 
the amplitude of the quasar-to-quasar errors listed in Table[3] 

Table|5]presents the master quasar templates used to estimate 
the cross-correlation redshifts. Col. 1 lists the rest-frame wave- 
length (A). Cols. 2, and 3 include the relative flux (per unit wave- 
length) and the number of spectra contributing for the master tem- 
plate, while Cols. 4 and 5 provide the same information for the 
FD-quasar template. The FD-quasar template does not extend quite 
as far to the blue and the flux column contains entries of '-999. 0' 
for wavelengths A< 1296.2 A. 

While the spectra should prove of use in the context of red- 
shift estimation, the templates are not suitable for studies of quasar 
spectral energy distributions, where care must be taken in defining 
the large-scale shape of such composite spectra. 
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2 3 

HW Redshift 

Figure 13. Redshift differences, A^/(l + z), between tlie final HW- 
redsliifts and Princeton redsliifts, as a function of redshift. Note the large 
range on the y-axis. Data for 90979 quasars are included (1306 lie outside 
the y-axis range plotted). 



9 DISCUSSION 

The approach adopted in this paper to the question of deriving red- 
shifts with a common zero-point over an extended dynamic range 
in redshift, and hence involving disjoint spectral wavelength cover- 
age, differs from that normally employed. The majority of studies 
to date have focussed on the parameterisation of the rest-frame cen- 
troid differences between the strongest emission lines present in the 
quasar spectra (e.g. Appendix A of Shen et al. 2007. for a recent 
example). Use of the cross-correlation redshifts directly, bypasses 
many of the difficulties associated in providing reliable, repro- 
ducible, parameterisations of low SNR, asymmetric, often blended, 
emission lines present on 'continua' that also show significant vari- 
ation from quasar to quasar. The resultant quasar-to-quasar disper- 
sion and the internal reproducibility of the new HW-redshifts rep- 
resent significant improvements over even the most careful studies 
utilising individual emission features. 

Systematic, luminosity-dependent relative emission line shifts 
have not featured in many previous studies of the quasar popula- 
tion. In part, the lack of such work may reflect the difficulty of 
performing such studies prior to the availability of the mo re recent 
SDSS Data Releases. An exception is the work of Rich ards et al.l 
( I2OO2I) who find a clear relationship between emission line cen- 
troid shifts, of ex actly the type discuss ed here, and emission line 
equivalent width. [Richards et al.l ( l2002h note that the line equiva- 
lent width is directly related to quasar absolute magnitude via the 
Baldwin Effect cBaldwin.1977,) . 



cedure/algorithm employed. However, comparison of the new HW- 
redshifts with the Princeton determinations (Fig. I13t still shows 
large (~600kms^^at 2:~2.2) systematic differences. 

The two evident 'jumps' in the relationship between the HW- 
and SDSS- estimates occur as cross-correlation redshifts includ- 
ing the Mg II A2799 emission line become important (at z~0.8) 
and where Mg II moves beyond the red limit of the SDSS spec- 
tra (at 2;~2. 1). The behaviour can be traced directly to differ- 
ences in the new quasar template and that of Vand en Berk et al.l 
|200ll)- Fig.nij shows the excellent agreement between the com- 
posites at optical wavelengths where emission lines, including H/3 
and [O III] AA4960,5008, dominate the redshift determinations (ei- 
ther directly, via emission line locations, or, through the contri- 
bution of emission lines to the cross-correlation signal). In con- 
trast. Fig. 114b . illustrates the significant difference in the location 
of the Mg II A2799 emission line in the two composites. At the ac- 
curacy levels of interest, absolute wavelength 'centroids' of broad 
emission lines in quasar spectra are dominated by the particular 
scheme used to define the associated 'continuum' and the height 
above the continuum used to define the 'line'. However, the cen- 
troid of the portion of the Mg II emission line above half the peak 
height is 1.2±0. lA blue r in the HW-template compared to the 
IVanden Berk et alj I2OO1I) template. The line centroid moves red- 
ward as increasingly large fractions of the line wings are included 
but at the half peak height level the HW-template line centroid is 
only 0.4±0.1 A (~45 km s~^) redward of the rest-frame reference 
value of 2798.75 A, derived from the Mg II components in the ratio 
2:1. 

The strongest 'jump' in the relation between the HW- and 
Princeton-redshifts in Fig. [T3] at z~2.2, derives fundamentally 
from the large systematic trends in the ratio of Mgll to 'Cm] ' 
emission line locations as a function of quasar absolute magni- 
tude (Fig. |5](. The origin of the effect is primarily the change in 
the ratio of C III] A 1908 to Si III] A 1892 (see Fig.fHlfor illustration 
in the context of FIRST-detected quasars) . A di rect comparison of 
the HW-template and Vanden Berk et^ j200lh template is some- 
what misleading because the HW-redshifts are derived only fol- 
lowing the significant absolute magnitude dependent corrections. 
However, using any sensible definition of the emission line, the 
C III] -I- Si III] blend is significantly bluer in the HW-template than 
in the Vanden Berk et al. (2001) composite, producing the increase 
in the HW-redshifts at 2:~2.2. The Princeton redshifts then become 
progressively closer to the HW-redshifts as the C IV emission line 
(with its well-established increasing blue asymmetry at increas- 
ing quasar luminosity) dominates the Princeton determinations at 
higher redshifts. Recall though, that the C IV emission region does 
not contribute to the HW-redshifts, except in a very small percent- 
age of quasars. 



9.1 Comparison with Prince ton redshifts and the 
IVanden Berk et alj j200ll) quasar template 

Operationally, it is found that straightforward systematic correc- 
tions to quasar redshift estimates, as a function of quasar ab- 
solute magnitude, reduce the systematic trends as a function 
of redshift (absolute magnitude) to <30kms~^per unit redshift 
(<10kms~^per magnitude). Internal reproducibility represents a 
factor >2 improvement over the SDSS redshift determinations. 
The results presented above, combined with the form of the dif- 
ferences between the Princeton and SDSS redshifts (Section 
show that the origin of a significant proportion of the improve- 
ments achieved are due to differences in the cross-correlation pro- 



9.2 Associated C iv and Mg 11 absorbers as quasar redshift 
diagnostics 

The availability of the large SDSS quasar catalogues have stim- 
ulated new investigations into the physical origin of associated 
absorbers, particularly those evident through the presence of 
Civ and Mgll absorption (e.g.lNestor. Hamaim. & Hidalgoll200g: 
IWild et alj200^:IVanden Berk et alj20 08TX pre-requisite for such 
investigations is an estimate of the systemic quasar redshifts. Given 
the large intrinsic variation in the properties of the C IV emission 
line and the relative invariance of the Mg II emission line centroid. 
redshifts based on the location of the Mg II emission line are often 
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Figure 14. Rest-frame spectra of the new quasar composite compared to that of IVanden Ber k et"al Panel (a) shows the prominent emission line of 

Mgl l AA2796.35,2803.53, wit h the new quasar composite (solid line) and the Vanden Berk et aP <200ll) composite (dashed line). The lower plot shows the 
ratio IVanden Berk et alj j200l Vnew. The vertical dotted line indicates the wavelength 2798.75 A, derived from the Mgll components in the ratio 2:1. The 
significantly bluer location of Mg II emission in the new composite is evident. Panel (b) shows the same information for the rest-frame wavelength region 
including H/3 and [O III] A A4960,5008. The vertical dotted lines indicate the rest-frame wavelengths of 4862.68 4960.30 and 5008.24 A for H,3 and [O III] . 
While the new composite possesses slightly stronger [O III] emission there is no evidence for any detectable offset in the emission line locations of H/3 or 
[OIII]. 
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Figure 15. Observed frequency distribution of redshift differences, /3=ti/c, 
for ~23 800 C IV absorbers using both SDSS- (blue) and HW-redshifts 
(red) for quasars with redshifts 1.55<2<3.5. The HW-redshift-based his- 
togram shows a much higher, better defined, peak, centred close to /3~0 
and a greatly reduced population of absorbers with positive fi values (i.e. 
^aba>^qao)- The centroid of the /3~0 component for the HW-redshifts 
shows no detectable shift over the entire redshift range of the quasars, 
1.55<2<3.5. 
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Figure 16. Observed frequency distribution of redshift differences, fi=v / c, 
for ~8750 Mgll absorbers using both SDSS- (blue) and HW-redshifts 
(red) for quasars with redshifts 0.45<2<2.1. The HW-redshift-based his- 
togram shows a higher, better defined, peak, centred close to /3~0 and a 
significantly reduced population of absorbers with positive /3 values (i.e. 
^aha^^qso). The centroid of the /3~0 component for the HW-redshifts 
shows no detectable shift over the entire redshift range of the quasars, 
0.45<z<2.1. 



employed in studies of both associated C IV and Mg II absorbers in 
quasars with redshifts 2^2.1. 

Both Mgll and CiV absorber catalogues are available from 
our investigation of absorber populations in SDSS quasars (e.g. 
IWild, Hewett. & Pettiniliooeh . Strong narrow absorbers are flagged 
and 'removed' from the quasar spectra prior to the calculation of 
the cross-correlation redshift determinations (Section |4l(. The new 
HW-quasar redshifts are thus essentially independent of the pres- 
ence of individual absorbers and a comparison of associated ab- 
sorbers velocity distributions, using both SDSS- and HW-redshifts, 
provides a powerful test of the redshift accuracy in an astrophysical 
context of considerable current interest. Fig.llSlshows the observed 



distribution of redshift differences EE P=v/c, for ~23 800 C IV ab- 
sorbers using both SDSS- and HW-redshifts for quasars with red- 
shifts 1.55<z<3.5. The differences in the distributions are strik- 
ing, with the HW-redshift-based histogram showing a much higher 
peak at /3~0 and a greatly reduced population of positive /3 val- 
ues (i.e. Zaba>Zqso)- The centroid of the /3~0 component shows 
no detectable shift over the entire redshift range of the quasars, 
1.55<2<3.5. 



The observed distributions are shown. No attempt has been made to 
calculate an absorber density by incorporating the redshift-path accessible 
as a function of /3. 
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Fig. [16] shows the equivalent distribution of redshift differ- 
ences, for ~8750 Mgll absorbers using both SDSS- and HW- 
redshifts for quasars with redshifts 0.45<2<2.1. The differences 
between the SDSS- and HW-redshifts at z<2.l are significantly 
smaller than for the quasars included in the C IV absorber sample. 
However, similar behaviour is evident to that seen in the C IV ab- 
sorbers, the /3~0 peak is significantly better defined using the HW- 
redshifts and the population of redshifted absorbers with positive /3 
greatly reduced. A more detailed consideration of the distribution 
of the Mg II absorber redshifts leads to an improved parameteri sa- 
tion of the various constituent absorber populations jWildl2009l) . 

To summarise, the distributions of the C IV and Mg II associ- 
ated absorber populations provide independent confirmation of the 
veracity of the new HW-redshifts. While beyond the scope of the 
present paper, investigations of associated absorber populations are 
in hand using the full SDSS DR7 Legacy Release spectroscopic 
database. 



10 CONCLUSIONS 

A systematic investigation of the relationship between differ- 
ent redshift estimation schemes for more than 91 000 quasars in 
the Sloan Digital Sky Survey (SDSS) Data Release 6 (DR6) is 
presented. Empirical relationships between redshifts based on i) 
Ca II H & K host galaxy absorption, ii) quasar [O II] A3728, iii) 
[Qui] AA4960,5008 emission, and iv) cross-correlation (with a 
master quasar template) that includes, at increasing quasar redshift, 
the prominent Mg II A2799, C III] A 1908 and C IV A 1549 emission 
lines, are established as a function of quasar redshift and luminos- 
ity. New redshifts in the resulting catalogue possess systematic bi- 
ases a factor of ~20 lower compared to the SDSS redshift values; 
systematic effects are reduced to the level of Az/(1 + 2)^10^* 
(30 km s^^) per unit redshift, or ^2.5 x 10^^ per unit absolute mag- 
nitude. 

It is important to realise that there will be systematic redshift 
trends present as a function of the quasar SEDs and the specific 
example of FIRST-detected quasars (Section[7]l provides an exam- 
ple, related to the radio-properties of the quasar SEDs. One of the 
primary motivations of this work is to facilitate further studies of 
SED-dependent systematic emission line properties, working from 
redshift estimates whose properties as a function of redshift and 
absolute magnitude are well understood. 

Equally important as the new redshift determinations, well- 
determined empirical estimates of the quasar-to-quasar dispersion 
in redshifts are available for each method of redshift estimation and 
a combined internal-l-population uncertainty is provided for every 
quasar in the catalogue. 

The improved redshifts and their associated errors have wide 
applicability in areas such as quasar absorption outflows, quasar 
clustering, quasar-galaxy clustering and proximity-effect determi- 
nations. 
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